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Hierarchical organization in complex networks 
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Many real networks in nature and society share two generic properties: they are scale- free and 
they display a high degree of clustering. We show that these two features are the consequence of 
a hierarchical organization, implying that small groups of nodes organize in a hierarchical manner 
into increasingly large groups, while maintaining a scale-free topology. In hierarchical networks the 
degree of clustering characterizing the different groups follows a strict scaling law, which can be 
used to identify the presence of a hierarchical organization in real networks. We find that several 
real networks, such as the World Wide Web, actor network, the Internet at the domain level and 
the semantic web obey this scaling law, indicating that hierarchy is a fundamental characteristic of 
many complex systems. 
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In the last few years an array of discoveries have re- 
defined our understanding of complex networks (for re- 
views see |^]). The availability of detailed maps, cap- 
turing the topology of such diverse systems as the cell 
[ 3 , [|, [|, |(| , the world wide web ]7j| , or the sexual network 
have offered scientists for the first time the chance to 
address in quantitative terms the generic features of real 
networks. As a result, we learned that networks are far 
from being random, but are governed by strict organiz- 
ing principles, that generate systematic and measurable 
deviations from the topology predicted by the random 
graph theory of Erdos and Renyi || [l0| , the basic model 
used to describe complex webs in the past four decades. 

Two properties of real networks have generated con- 
siderable attention. First, measurements indicate that 
most networks display a high degree of clustering. Defin- 
ing the clustering coefficient for node i with ki links as 
Cj = 2m/ki(ki — 1), where n, is the number of links 
between the ki neighbors of i, empirical results indicate 
that Ci averaged over all nodes is significantly higher for 
most real networks than for a random network of similar 
size jl], H[ pT[ . Furthermore, the clustering coefficient of 
real networks is to a high degree independent of the num- 
ber of nodes in the network (see Fig. 9 in (lj). At the 
same time, many networks of scientific or technological 
interest, ranging from the World Wide Web |t]] to biolog- 
ical networks j|, |J, ||, |(| have been found to be scale-free 
p2| , 13] ], which means that the probability that a ran- 
domly selected node has k links (i.e. degree k) follows 
P{k) ~ fc~ 7 , where 7 is the degree exponent. 

The scale-free property and clustering are not exclu- 
sive: for a large number of real networks, including 
metabolic networks |j, [I| , the protein interaction network 
[^|, |f|, the world wide web and even some social net- 
works 1 14, [l5|, |l6) the scale- free topology and high clus- 
tering coexist. Yet, most models proposed to describe 
the topology of complex networks have difficulty captur- 
ing simultaneously these two features. For example, the 
random network model ^, cannot account neither for 
the scale-free, nor for the clustered nature of real net- 



works, as it predicts an exponential degree distribution, 
and the average clustering coefficient, C(N), decreases 
as A^ 1 with the number of nodes in the network. Scale- 
free networks, capturing the power law degree distribu- 
tion, predict a much larger clustering coefficient than a 
random network. Indeed, numerical simulations indicate 
that for one of the simplest models 12, 13 the 



aver- 
age clustering coefficient depends on the system size as 
C(N) ~ A"" ' 75 U, g, significantly larger for large A 
than the random network prediction C(N) ~ N^ 1 . Yet, 
this prediction still disagrees with the finding that for 
several real systems C is independent of A" Q . 

Here we show that the fundamental discrepancy be- 
tween models and empirical measurements is rooted in a 
previously disregarded, yet generic feature of many real 
networks: their hierarchical topology. Indeed, many net- 
works are fundamentally modular: one can easily identify 
groups of nodes that are highly interconnected with each 
other, but have only a few or no links to nodes outside of 
the group to which they belong to. In society such mod- 
ules represent groups of friends or coworkers 17 ; in the 
WWW denote communities with shared interests jig, [l9| ; 
in the actor network they characterize specific genres or 
simply individual movies. Some groups are small and 
tightly linked, others are larger and somewhat less in- 
terconnected. This clearly identifiable modular organi- 
zation is at the origin of the high clustering coefficient 
seen in many real networks. Yet, models reproducing 
the scale-free property of real networks |jj || distinguish 
nodes based only on their degree, and are blind to node 
characteristics that could lead to a modular topology. 

In order to bring modularity, the high degree of clus- 
tering and the scale-free topology under a single roof, we 
need to assume that modules combine into each other in 
a hierarchical manner, generating what we call a hier- 
archical network. The presence of a hierarchy and the 
scale-free property impose strict restrictions on the num- 
ber and the degree of cohesiveness of the different groups 
present in a network, which can be captured in a quanti- 
tative manner using a scaling law, describing the depen- 
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dence of the clustering coefficient on the node degree. We 
use this scaling law to identify the presence of a hierarchi- 
cal architecture in several real networks, and the absence 
of such hierarchy in geographically organized webs. 



I. HIERARCHICAL NETWORK MODEL 

We start by constructing a hierarchical network model, 
that combines the scale-free property with a high degree 
of clustering. Our starting point is a small cluster of five 
densely linked nodes (Fig. j]a). Next we generate four 
replicas of this hypothetical module and connect the four 
external nodes of the replicated clusters to the central 
node of the old cluster, obtaining a large 25-node module 
(Fig. ^d) . Subsequently, we again generate four replicas 
of this 25-node module, and connect the 16 peripheral 
nodes to the central node of the old module (Fig. |l|c) , 
obtaining a new module of 125 nodes. These replication 
and connection steps can be repeated indefinitely, in each 
step increasing the number of nodes in the system by a 
factor five. 



(a) n=0, N=5 





(b) n=l, N=25 



(c) n=2, N=125 



FIG. 1: The iterative construction leading to a hierarchical 
network. Starting from a fully connected cluster of five nodes 
shown in (a) (note that the diagonal nodes are also connected 
- links not visible), we create four identical replicas, connect- 
ing the peripheral nodes of each cluster to the central node of 
the original cluster, obtaining a network of N — 25 nodes (b). 
In the next step we create four replicas of the obtained cluster, 
and connect the peripheral nodes again, as shown in (c), to 
the central node of the original module, obtaining a, N — 125 
node network. This process can be continued indefinitely. 

Precursors to the model described in Fig.jl] have been 
proposed in Ref. |^0| and extended and discussed in Ref. 
pll |22| as a method of generating deterministic scale- free 
networks. Yet, it was believed that aside from their deter- 
ministic structure, their statistical properties are equiv- 
alent with the stochastic models that are often used to 
generate scale-free networks. In the following we argue 
that such hierarchical construction generates an architec- 



ture that is significantly different from the networks gen- 
erated by traditional scale- free models. Most important, 
we show that the new feature of the model, its hierarchi- 
cal character, are shared by a significant number of real 
networks. 

First we note that the hierarchical network model 
seamlessly integrates a scale- free topology with an inher- 
ent modular structure. Indeed, the generated network 
has a power law degree distribution with degree expo- 
nent 7=1 + In5/ln4 = 2.161 (Fig.||a). Furthermore, 
numerical simulations indicate that the clustering coeffi- 
cient, C ~ 0.743, is independent of the size of the net- 
work (Fig. ||c;). Therefore, the high degree of clustering 
and the scale-free property are simultaneously present in 
this network. 

The most important feature of the network model of 
Fig. [j], not shared by either the scale-free (lj, [ll| or ran- 
dom network models [To), is its hierarchical architec- 
ture. The network is made of numerous small, highly 
integrated five node modules (Fig.fjla), which are assem- 
bled into larger 25-node modules (Fig. |l|b) . These 25- 
node modules are less integrated but each of them is 
clearly separated from the other 25-node modules when 
we combine them into the even larger 1 25-node modules 
(Fig. [ljs). These 125-node modules are even less cohesive, 
but again will appear separable from their replicas if the 
network expands further. 

This intrinsic hierarchy can be characterized in a quan- 
titative manner using the recent finding of Dorogovtsev, 
Goltsev and Mendes (2l| that in deterministic scale-free 
networks the clustering coefficient of a node with k links 
follows the scaling law 



C(k) ~ k- 



(1) 



We argue that this scaling law quantifies the coexis- 
tence of a hierarchy of nodes with different degrees of 
clustering, and applies to the model of Fig. 0a-c as well. 
Indeed, the nodes at the center of the numerous 5-node 
modules have a clustering coefficient C = 1. Those at the 
center of a 25-node module have k — 20 and C = 3/19, 
while those at the center of the 125-node modules have 
k = 84 and C — 3/83, indicating that the higher a node's 
degree the smaller is its clustering coefficient, asymptot- 
ically following the l/k law (Fig.|2|b). In contrast, for 
the scale-free model proposed in Ref. |l^| the clustering 
coefficient is independent of k, i.e. the scaling law (|l|) 
does not apply (Fig.||b). The same is true for the ran- 
dom H |l0| or the various small world models Eg) , 
for which the clustering coefficient is independent of the 
nodes' degree. 

Therefore, the discrete model of Fig. [I] combines within 
a single framework the two key properties of real net- 
works: their scale-free topology and high modularity, 
which results in a system-size independent clustering co- 
efficient. Yet, the hierarchical modularity of the model 
results in the scaling law ([!]), which is not shared by 
the traditional network models. The question is, could 
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FIG. 2: Scaling properties of the hierarchical model shown in Fig. |lj (N = 5 7 ). (a) The numerically determined degree 
distribution. The assymptotic scaling, with slope 7 = 1 + In 5/ hi 4, is shown as a dashed line, (b) The C(k) curve for the 
model, demonstrating that it follows Eq. (|l|). The open circles show C{k) for a scale-free model jijj of the same size, illustrating 
that it does not have a hierarchical architecture, (c) The dependence of the clustering coefficient, C, on the size of the network 
N. While for the hierarchical model C is independent of N (♦), for the scale-free model C(N) decreases rapidly (o). 



hierarchical modularity, as captured by this model, char- 
acterize real networks as well? 



II. HIERARCHICAL ORGANIZATION IN REAL 
NETWORKS 

To investigate if such hierarchical organization is 
present in real networks we measured the C(k) function 
for several networks for which large topological maps are 
available. Next we discuss each of these systems sepa- 
rately. 

Actor Network: Starting from the www . IMDB . com 
database, we connect any two actors in Hollywood if they 
acted in the same movie, obtaining a network of 392,340 
nodes and 15,345,957 links. Earlier studies indicate that 
this network is scale-free with an exponential cutoff in 
P(k) for high k §J|, ||[ As Fig.|a indicates, we 

find that C(k) scales as indicating that the net- 

work has a hierarchical topology. Indeed, the majority 
of actors with a few links (small k) appear only in one 
movie. Each such actor has a clustering coefficient equal 
to one, as all actors the actor has links to are part of 
the same cast, and are therefore connected to each other. 
The high k nodes include many actors that acted in sev- 
eral movies, and thus their neighbors are not necessarily 
linked to each other, resulting in a smaller C(k). At high 
k the C(k) curve splits into two branches, one of which 
continues to follow Eq. ([]]), while the other saturates. 
One explanation of this split is the decreasing amount of 
datapoints available in this region. Indeed, in the high k 
region the number of nodes having the same k is rather 
small. If one of these nodes corresponds to an actor that 
played only in a few movies with hundreds in the cast, it 
will have both high k and high C, considerably increasing 
the average value of C(k). The k values for which such a 



high C nodes are absent continue to follow the curve, 
resulting in jumps between the high and small C values 
for large k. For small k these anomalies are averaged out. 

Language network: Recently a series of empirical re- 
sults have shown that the language, viewed as a network 
of words, has a scale- free topology p(l 12?], 28, p9). Here 



we study the network generated connecting two words to 
each other if they appear as synonyms in the Merriam 
Webster dictionary [^?J . The obtained semantic web has 
182,853 nodes and 317,658 links and it is scale-free with 
degree exponent 7 = 3.25. The C(k) curve for this lan- 
guage network is shown in Fig. [3)3, indicating that it fol- 
lows (Q), suggesting that the language has a hierarchical 
organization. 

World Wide Web: On the WWW two documents are 
connected to each other if there is an URL pointing from 
one document to the other one. The sample we study, ob- 
tained by mapping out the www . nd . edu domain JtJ , has 
325,729 nodes and 1,497,135 links, and it is scale-free 
with degree exponents 7 ou t = 2.45 and 7i n = 2.1, charac- 
terising the out and in-degree distribution, respectively. 
To measure the C{k) curve we made the network undi- 
rected. While the obtained C(k), shown in Fig.||c, does 
not follow as closely the scaling law (]]]) as observed in the 
previous two examples, there is clear evidence that C{k) 
decreases rapidly with k, supporting the coexistence of 
many highly interconnected small nodes with a few larger 
nodes, which have a much lower clustering coefficient. 

Indeed, the Web is full of groups of docu- 
ments that all link to each other. For example, 
www . nd . edu/~networks, our network research dedicated 
site, has a high clustering coefficient, as the documents 
it links to have links to each other. The site is one of 
the several network-oriented sites, some of which point 
to each other. Therefore, the network research commu- 
nity still forms a relatively cohesive group, albeit less in- 
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terconnected than the www.nd. edu/~networks site, thus 
having a smaller C . This network community is nested 
into the much larger community of documents devoted 
to statistical mechanics, that has an even smaller cluster- 
ing coefficient. Therefore, the fc-dependent C{k) reflects 
the hierarchical nesting of the different interest groups 
present on the Web. Note that C(k) ~ fc _1 for the 
WWW was observed and briefly noted in Ref. |30|. 




FIG. 3: The scaling of C{k) with k for four large networks: 
(a) Actor network, two actors being connected if they acted in 
the same movie according to the www . IMDB . com database, (b) 
The semantic web, connecting two English words if they are 
listed as synonyms in the Merriam Webster dictionary | pj . 
(c) The World Wide Web, based on the data collected in Rcf. 
M. (d) Internet at the Autonomous System level, each node 
representing a domain, connected if there is a communication 
link between them. The dashed line in each figure has slope 
-1, following Eq. (@). 

Internet at the AS level: The Internet is often studied 
at two different levels of resolution. At the router level we 
have a network of routers connected by various physical 
communication links. At the interdomain or autonomous 
system (AS) level each administrative domain, composed 
of potentially hundreds of routers, is represented by a sin- 
gle node. Two domains are connected if there is at least 
one router that connects them. Both the router and the 
domain level topology have been found to be scale-free 
pl[ . As Fig.^d shows, we find that at the domain level 
the Internet, consisting of 65,520 nodes and 24,412 links 
HU, has a hierarchical topology as C(k) is well approxi- 
mated with (0). The scaling of the clustering coefficient 
with k for the Internet was earlier noted by Vazquez, 
Pastor-Satorras and Vespignani (VPSV) |33, |I|, who 
observed C(k) ~ k~ 0J5 . VPSV interpreted this finding, 
together with the observation that the average nearest- 
neighbor connectivity also follows a power-law with the 
node's degree, as a natural consequence of the stub and 
transit domains, that partition the network in a hier- 
archical fashion into international connections, national 



backbones, regional networks and local area networks. 

Our measurements indicate, however, that some real 
networks lack a hierarchical architecture, and do not obey 
the scaling law (Q). In particular, we find that the power 
grid and the router level Internet topology have a k in- 
dependent C(k). 

Internet at the router level: The router level Inter- 
net has 260,657 nodes connected by 1,338,100 links [p5| . 
Measurements indicate that the network is scale-free 
lHI, I]! with degree exponent 7 = 2.23. Yet, the C{k) 
curve (Fig.|]a), apart from some fluctuations, is largely 
independent of k, in strong contrast with the C(k) ob- 
served for the Internet's domain level topology (Fig.||d), 
and in agreement with the results of VPSV ]33|, |4j, who 
also note the absence of a hierarchy in router level maps. 

Power Grid: The nodes of the power grid are genera- 
tors, transformers and substations and the links are high 
voltage transmission lines. The network studied by us 
represents the map of the Western United States, and 
has 4,941 nodes and 13,188 links pi] . The results again 
indicate that apart from fluctuations, C(k) is indepen- 
dent of k. 

It is quite remarkable that these two networks share a 
common feature: a geographic organization. The routers 
of the Internet and the nodes of the power grid have a 
well defined spatial location, and the link between them 
represent physical links. In contrast, for the examples 
discussed in Fig. || the physical location of the nodes was 
either undefined or irrelevant, and the length of the link 
was not of major importance. For the router level Inter- 
net and the power grid the further are two nodes from 
each other, the more expensive it is to connect them 
|p6[ . Therefore, in both systems the links are driven by 
cost considerations, generating a distance driven struc- 
ture, apparently excluding the emergence of a hierarchi- 
cal topology. In contrast, the domain level Internet is less 
distance driven, as many domains, such as the AT&T do- 
main, span the whole United States. 

In summary, we offered evidence that for four large 
networks C(k) is well approximated by C(k) ~ fc _1 , in 
contrast to the fc-independent C(k) predicted by both the 
scale-free and random networks. In addition, there is ev- 
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FIG. 4: The scaling of C(k) for two large, non-hierarchical 
networks: (a) Internet at router level [p5| . (b) The power 
grid of Western United States. The dashed line in each figure 
has slope —1, while the solid line corresponds to the average 
clustering coefficient. 
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idence for similar scaling in the metabolism [ B7j and pro- 
tein interaction networks 38 . This indicates that these 



networks have an inherently hierarchical organization. In 
contrast, hierarchy is absent in networks with strong geo- 
graphical constraints, as the limitation on the link length 
strongly constraints the network topology. 



III. STOCHASTIC MODEL AND 
UNIVERSALITY 

The hierarchical model described in Fig. |l| predicts 
C(k) ~ k _1 , which offers a rather good fit to three of 
the four C(k) curves shown in Fig.||. The question is, 
is this scaling law (0) universal, valid for all hierarchical 
networks, or could different scaling exponent characterize 
the scaling of C(k)l Defining the hierarchical exponent, 
/?, as 

k-P 



C{k) 



(2) 



is (3 = 1 a universal exponent, or it's value can be changed 
together with 7? In the following we demonstrate that 
the hierarchical exponent f3 can be tuned as we tune some 
of the network parameters. For this we propose a stochas- 
tic version of the model described in Fig. 

We start again with a small core of five nodes all con- 
nected to each other (Fig.[l]a) and in step one (n = 1) 
we make four copies of the five node module. Next, we 
randomly pick a p fraction of the newly added nodes and 
connect each of them independently to the nodes belong- 
ing to the central module. We use preferential attach- 
ment |l3| to decide to which central node the selected 
nodes link to. That is, we assume that the probability 
that a selected node will connect to a node i of the cen- 
tral module is ki / • kj , where fc; is the degree of node 
i and the sum goes over all nodes of the central module. 
In the second step (n = 2) we again create four identical 
copies of the 25-node structure obtained thus far, but we 
connect only a p 2 fraction of the newly added nodes to 
the central module. Subsequently, in each iteration n the 
central module of size 5™ is replicated four times, and in 
each new module a p n fraction will connect to the cur- 
rent central module, requiring the addition of (5p) n new 
links. 

As Fig. [3] shows, changing p alters the slope of both 
P(k) and C(k) on a log-log plot. In general, we find that 
increasing p decreases the exponents 7 and j3 (Fig. ^>,d). 
The exponent (3 = 1 is recovered for p = 1, i.e. when all 
nodes of a module gain a link. While the number of links 
added to the network changes at each iteration, for any 
p < 1 the average degree of the infinitely large network 
is finite. Indeed, the average degree follows 



(k), 



'3 
r, V2 



1-p 



n+l 



1-p 



(3) 



which is finite for any p < 1. 

Interestingly, the scaling of C(k) is not a unique prop- 
erty of the model discussed above. A version of the 
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FIG. 5: The scaling properties of the stochastic model, (a) 
The degree distribution for different p values, indicating that 
P(k) follows a power law with a p dependent slope, (b) The 
dependence of the degree exponent 7 on p, determined by 
fitting power laws to the curves shown in (a). The exponent 
7 appears to follow approximately j(p) ~ 1/p (dashed line), 
(c) The C(k) curve for different p values, indicating that the 
hierarchical exponent /3 depends on p. (d) The dependence 
of (5 on the parameter p. The simulations were performed for 
N = 5 7 (78,125) nodes. 



model, where we keep the fraction of selected nodes, p, 
constant from iteration to iteration, also generates p de- 
pendent (3 and 7 exponents. Furthermore, recently sev- 
eral results indicate that the scaling of C(k) is an intrin- 
sic feature of several existing growing networks models. 
Indeed, aiming to explain the potential origin of the scal- 
ing in C(k) observed for the Internet, VSPV note that 
39, |o) displays a C(k) that appears 
While there is no analytical evidence 
numerical results [B3l p54fl suggest 



the fitness model 
to scale with k. 

-0 



for C(k) ~ k~P yet, 
that the presence of fitness does generate a hierarchi- 
cal network architecture. In contrast, in a recent model 
proposed by Klemm and Eguiluz there is analytical evi- 
dence that the network obeys the scaling law (Q) Q . In 
their model in each time step a new node joins the net- 
work, connecting to all active nodes in the system. At 
the same time an active node is deactivated with prob- 
ability p ~ The insights offered by the hierarchical 
model can help understand the origin of the observed 
C(k) ~ fc _1 . By deactivating the less connected nodes a 
central core emerges to which all subsequent nodes tend 
to link to. New nodes have a large C and small k, thus 
they are rapidly deactivated, freezing into a large C state. 
The older, more connected, surviving nodes are in con- 
tact with a large number of nodes that have already dis- 
appeared from the active list, and they have small C [p2| . 

Finally, Szabo, Alava and Kertesz have developed a 
rate equation method to systematically calculate C(k) 
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for evolving networks models j43fl . Applying the method 
to a model proposed by Holme and Kim mM to enhance 
the degree of clustering coefficient C seen in the scale- 
free model [|l2|, they have shown that the scaling of C(k) 
depends on the parameter p, which governs the rate at 
which new nodes connect to the neighbors of selected 
nodes, bypassing preferential attachment. As for p = 
the Holme-Kim model reduces to the scale-free model, 
Szabo, Alava and Kertesz find that in this limit the scal- 
ing of C(k) vanishes. These models indicate that sev- 
eral microscopic mechanisms could generate a hierarchi- 
cal topology, just as several models are able to create a 
scale- free network [[l]. ^| . 



network. 
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IV. DISCUSSION AND OUTLOOK 



The identified hierarchical architecture offers a new 
perspective on the topology of complex networks. In- 
deed, the fact that many large networks are scale-free 
is now well established. It is also clear that most net- 
works have a modular topology, quantified by the high 
clustering coefficient they display. Such modules have 
been pro pos ed to be a fundamental feature of biological 
systems |37|, but have been discussed in the con- 



text of the WWW 1 18, |4q] , and social networks as well 
[f7j . The hierarchical topology offers a new avenue for 
bringing under a single roof these two concepts, giving a 
precise and quantitative meaning for the network's mod- 
ularity. It indicates that we should not think of modu- 
larity as the coexistence of relatively independent groups 
of nodes. Instead, we have many small clusters, that are 
densely interconnected. These combine to form larger, 
but less cohesive groups, which combine again to form 
even larger and even less interconnected clusters. This 
self-similar nesting of different groups or modules into 
each other forces a strict fine structure on real networks. 

Most interesting is, however, the fact that the hier- 
archical nature of these networks is well captured by a 
simple quantity, the C(k) curve, offering us a relatively 
straightforward method to identify the presence of hier- 
archy in real networks. The law (|l|) indicates that the 
number and the size of the groups of different cohesive- 
ness is not random, but follow rather strict scaling laws. 

The presence of such a hierarchical architecture rein- 
terprets the role of the hubs in complex networks. Hubs, 
the highly connected nodes at the tail of the power law 
degree distribution, are known to play a key role in keep- 
ing complex networks together, playing a crucial role 
from the robustness of the network M, H9| to the spread 
of viruses in scale- free networks |5^] . Our measurements 
indicate that the clustering coefficient characterizing the 
hubs decreases linearly with the degree. This implies that 
while the small nodes are part of highly cohesive, densely 
interlinked clusters, the hubs are not, as their neighbors 
have a small chance of linking to each other. Therefore, 
the hubs play the important role of bridging the many 
small communities of clusters into a single, integrated 
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